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Abstract 

We have found a simple criterion which ahows for the straightfor- 
ward determination of the order-disorder critical temperatures. The 
method reproduces exactly results known for the two dimensional 
Ising, Potts and Z{N < 5) models. It also works for the Ising model on 
the triangular lattice. For systems which are not selfdual our propo- 
sition remains an unproven conjecture. It predicts Pc = 0.2656... for 
the two coupled layers of Ising spins. Critical temperature of the three 
dimensional Ising model is related to the free energy of the two layer 
Ising system. 
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Classical method of determining critical temperature in statistical physics 
consists of locating singularities of the largest eigenvalue of the corresponding 
transfer matrix T. To this end one has to study the high power of T, e.g. 
with L being the linear size of the system. This amounts to investigate the d 
dimensional euclidean system in its full complexity. On the other hand it is 
conceivable that the information about the phase transition is also encoded 
in all other eigenvalues, i.e. that the whole spectrum of the transfer matrix is 
sensitive to its location. This observation is confirmed by the famous example 
of the two dimensional Ising model where all eigenvalues show the extremal 
behaviour at = /3c Therefore we propose to search for the extremum 
of some simpler function which characterizes the system, which however is 
not dominated solely by the largest eigenvalue. The advantage of such an 
approach is that one may infer a nontrivial information from " characteristic" 
functions which are much easier to calculate. To be specific we propose to 
study the following characteristic function of a d dimensional system 

and analogous higher, moments of the transfer matrix. Following simple 
properties of p can be easily proven. 

a) p(0) = 1, 

h) p(oo) = l, (2) 
c) p(/3)>l. 

To see a) consider the normalized second moment at finite L 

rM = -jr;^ = (3) 

where Zi and Z2 are the partition functions of the d — 1 dimensional system 
and of the two coupled d — 1 dimensional systems respectively. For simplic- 
ity we shall use the terminology of the d=2 spin system. Hence Zi and Z2 
describe one dimensional spin chain and two coupled chains of spins. Peri- 
odic boundary conditions are implied in both directions, even for the single 
chain. With these definitions a) follows immediately from the observation 
that Z2{P = 0) = Af2 = (A/i)^, where A/i,A/'2 denote the total number of 
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microscopic states of single and double chain. At low temperatures {P 
fully ordered states dominate, hence 



oo 



lim g'' 



1. 



oo) 



(4) 



Where g denotes degeneracy factor of the ordered states. The last equality 
requires finite g (and d > 1), therefore b) holds only for systems with the 
discrete internal symmetry. Finally the property c) follows directly from the 
positivity of the eigenvalues of the transfer matrix since 



> 1. 



(5) 



Characteristic function, Eq.(|l|), is sensitive to the whole spectrum of the 
transfer matrix. However in view of our earlier discussion supplemented with 
the properties a) — c) it is natural to expect that the maximum of p occurs 
at the phase transition point. 
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Surprisingly this simple proposition is true in many, sometimes nontrivial, 
cases. We shall discuss them in the order of increasing complexity. 

Two dimensional Ising model. For d = 2, calculation of p{(3) is easily 
reduced to solving a straightforward two spin problem. Indeed the partition 
functions Zi[i = 1,2) are readily written as 



(7) 



where are the transfer matrices propagating one/two spins horizontally, 
c.f. Fig.l. 



Explicitly for Ising spins 
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We have chosen the interaction energy 

m = - E 



(9) 



<nm> 
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and s„ = ±1. In the thermodynamical limit the largest eigenvalues of Tj 
dominate and we get 

PIs^n,2m = = X\X + 1)' 

t2max i(x2 + 1)2 + ^ \{x^ - 1)4 + x^{x^ + 1)2 

Since the L ^ oo limit was already performed, p(/3), as given by Eg. ([T0|) , is 
the nontrivial characteristic of the infinite system. As conjectured p(/5) has 
a single maximum at Xc = 1 + which corresponds to the famous Onsager 
value ilQ. 

Potts and equivalent Z{N) models. The energy of the q state Potts model 
is given by Eq.(^) with s„ assuming q different values. Transfer matrices for 
one and two spin system are again simple 

<s|Ti|s'> = expp{Ss,s' + l), (11) 

< Si, S2\T2\s[, S'2 > = expP{6s^^s2 + + Ss2,s'2 + ^s\,sO- 

For q = 3 the diagonalization is tractable and the final result for the 
characteristic function reads 



PPottsS 



2xHx + 2) 



■4 + 3x2 + 22; + 3 + ^a;8 + 2x^ - 4x5 + 27x^ + 28x^ + 6x2 ^ + 9 



X' 



Again it has a single maximum at Xc = 1 + which agrees with the known 
location of the transition temperature 0. 

For g = 4, 5 we have used standard numerical methods to diagonalize 
transfer matrices. In both cases p(x) has a sigle maximum located at Xc = 
1 + y/g in accord with known results. For g > 4 the transition is first order 
[|] and consequently our method seems to apply to both kinds of transitions. 
We have not attempted algebraic diagonalization of Tj for arbitrary q. Low 

(A^ < 5) Z[N) (clock) models are equivalent to Potts3 (N=3) and Ising 
(N=4) systems. Not surprisingly the maximum of p{P) again agrees with 
known results. We obtain xf^^^ = (1 + VS)^ and xf*^*^) = 1 + V2 within the 
accuracy of our numerical procedures. 

^Note the difference by a factor of 2 which is caused by our choice of the Potts-hke 
interactions in Eq.(^. 

^AU algebraic calculations were done with the aid of MATHEMATICA. 
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All previously considered systems were self dual. Therefore one may jus- 
tifiably wonder if our principle is not yet another manifestation of the selfd- 
uality. Next example demonatrates that the "maximum rule", Eq.(^), is at 
least more general than the simple duality. 

Ising model on a triangular lattice. Transfer matrices Ti 2 (see Fig. 2.) 

<si|Ti|s'i> = exp/?(2sis'i + 1), (12) 

< Sl, S2\T2\s[, 4 > = expP{sis'i + S2S2 + S1S2 + S[S2 + S1S2 + S2s'i), 

can be simply diagonalized. We get 

tr" = -^(3 + a;^ + - 2x8 + 16x^ + 1). (14) 

The maximum of the characteristic function is located at Xmax = 3^. This 
agrees with the critical temperature first derived for this system by On- 
sager 0. One should remember however that Ising models on triangular and 
hexagonal lattices are interrelated via the duality and star-trangle relations, 
hence effectively there exists a symmetry which determines transition points 
in both systems ||^. Our maximum rule could in principle be a consequence 
of such a symmetry in this case. Next application provides more stringent 
test of this possibility. 

Two layer Ising model. This system consists of the two planes of Ising 
spins coupled by the nearest neighbour ferromagnetic interaction, also along 
the vertical (between the planes) direction. Periodic boundary conditions 
are assumed in all three directions, which amounts to doubling the strengh 
of interaction between the planes. The system in not selfdual and to our 
knowledge no other more complicated symmetries are known. Consequently 
its critical temperature was never derived. On the contrary our method pro- 
vides relatively simple analytic predictions for f3max- We proceed analogously 
to the previous cases. Transfer matrices Ti 2 propagate now states of two and 
four spins respectively (see Fig. 3). Ti reads 

< siS2\Ti\s\s'2>= ey^pP{siS2 + sis[ + 3282 + s[s'2 + 2), Sj = ±1, (15) 

■^From now on we use the standard Ising energy E{s) = — X]<i7> ■SjSj. 
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and its largest eigenvalue is 



2 

^rnax _ ^ (3.4 ^ 3 + + Vx^ + 14 + X^^) . (16) 

Matrix elements of are 

< ■SiS2S3S4|T2|SiS24'54 >= exp/5(SiS2 + S2S3 + S3S4 + S4S1) 

exp i5{sis\ + S2S2 + S3S3 + S4S4) exp i5{s\s'2 + 44 + -s's^ + "Sl-s'i)- 

Diagonalisation of this 16 x 16 matrix is simplified noting that conserves 
the [/-parity, [72,?/] = 0, [/ = Vu=\ '^f- The largest eigenvalue belongs to the 
U = +1 sector. Final expression is little more complicated 



= ^^i^^3 + ^v^+:^V^2 + (l + a:S)g3y^, (17) 



4x12 4x12 2V2. 



X 



q^(x) = x^2- 4x2^ + 70x^^-4x^ + 1, (18) 
g2(x) = x^° - 2x^^ + 5x^2 + 26x^^ + 4x^° + 26x^*^ + 5x^-2x^ + 1, 
q^(x) = x^^ - 2x^2 + 6x^ -2x^ + 1. 

Resulting characteristic function p(/3) is shown in Fig.4. It has the single 
maximum located at 

Prnax = 0.2656... . (19) 

According to our proposition, Eq.(^), this gives the transition temperature 
of the two layer Ising system. For comparison: /?^'^*"^2 ~ 0.4407... and 
pising3 0.221652(3) 0. Confronting this number with the results from 
MC simulations would provide the crucial test of our hypothesis. It would 
be also very interesting to search for the "generalized self duality" - the in- 
variance which would assure existence of a single maximum at the transition 
point. An attractive possibility is to use the characteristic function given by 
Eqs.(|16],|18|) to define such a mapping. This assumption has many verifi- 



able consequences. For example, all higher normalized moments of T should 
respect the same symmetry. 

Three dimensional Ising model. According to our proposition the tran- 
sition temperature of the d dimensional system is determined by the (3 de- 
pendence of the free energies 2 = — log Zi 2 of the corresponding d — 1 
dimensional systems. In particular, the problem of finding /^^ for the three 
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3 
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6 7 




0.3317 


0.3067 


0.2938 


0.2859 0.2698 



Table 1: Table I. Volume dependence of the pseudocritical temperature for 
the three dimensional Ising model. 



dimensional Ising model would be reduced to finding the free energy of the 
two coupled layers of the Ising spins. Indeed in this case the full transfer ma- 
trix T propagates the whole plane of spins, say, vertically, while the reduced 
transfer matrices Ti 2 propagate one (two) rows horizontally. No analytic 
solution of this system exists up to date Q However exact expressions for 
the complete partition functions Zi^2{P,L) in the finite volume are available 
for not so large L |^. We have therefore calculated the exact locations 
PmaxiL) of the ratio ri{(3), c.f. Eq. (|), which define the pseudocritical tem- 
perature at finite volume, c.f. Table I. 

In the thermodynamical limit f3max{L) should converge to the true tran- 
sition temperature. Even though the available range of L values is rather 
limited one sees the proper trend in the L dependence. Our values definitely 
move toward f3^-''^^9^ which was quoted above. One needs larger sizes to test 
quantitatively the L dependence against the finite size scaling predictions. 

Variety of approximate methods (Monte Carlo, high temperature expan- 
sion) can be also employed to test predictions of Eq.(^, in this case. 

Limits of applicability. We have also investigated situations where the 
maximum rule does not work. The regularity emerging from this study indi- 
cates that the method does not apply to systems with more than two phases. 
We have calculated the characteristic function for the variety of models with 
the intermediate Kosterlitz-Thouless (KT) phase. The maximum always 
occured inside the KT region. This phenomenon was found for Z{N) mod- 
els with N=5-19, and for the icosahedron and dodecahedron models in two 
dimensions. For the 0(2) model Pmax = 1.35... well above the known MC esti- 
mate for the transition between the disordered and KT phases (3c = 1.1197(5) 
[p!0| , |lT| . Interpretation of Pmax located inside the KT phase remains an open 



and interesting problem. On the other hand in the two dimensional 0(3) 
model p{(3) does not have any maximum in accord with the common wisdom 

^In our previous application we have derived the transition temperature only. 
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about the lack of a phase transition in this model. Note that for 0(2) and 
0(3) models the property c) is not satisfied and indeed calculated p{oo) > 1. 
Nevertheless the characteristic function has a maximum for 0(2) while it is 
monotonic for 0(3). This parallels the difference of the phase structures of 
both models. 

Summary. We have found a surprisingly simple criterion for locating 
order-disorder transition. It is exact for selfdual models. The method allows 
for the analytic calculation in variety of more complicated systems. In par- 
ticular we give the analytic estimate of the critical temperature of the two 
layer Ising system. Monte Carlo check of this prediction should be the first 
step towards more advanced applications. 

Our proposition reduces determination of the critical temperature of the 
three dimensional Ising model to finding the (3 dependence of the free energy 
of the two coupled planes of Ising spins. 

While we are lacking the complete proof of our hypothesis in general 
case, many approximated methods can be employed to test it in specific 
applications. High temperature expansion is one interesting possibility. 

One can reformulate the maximum principle in other equivalent ways. 
Differentiating the logarithm of Eq.(^) gives as the condition for the maxi- 
mum 

U2i/3c) = uii/3,), (20) 

where Mi 2 denotes the density of the internal energy for the one (two) layer 
system. Analogous relations follow from applying our proposition to higher 
moments of transfer matrix. They all say that at the bulk {d = 3 say) critical 
temperature internal energies of the interacting and noninteracting planes 
are equal. According to our earlier discussion this statement follows from 
duality for the two dimensional Ising and Potts models with planes replaced 
by the spin chains. Perhaps the most interesting formulation of Eq. (|20|) for 
higher moments of the transfer matrix results in the limit of infinite number 
of planes. Then we recover the original two dimensional system and our 
criterion reads 

uM'^) = ^<i-M''^), (21) 

where denotes the internal energy density of the d dimensional system, 
and Pj:"^^ corresponds to its crtitical temperature. As emphasized before, this 
statement follows from self duality in the case of the two dimensional Ising 
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model. It can be also checked directly. Indeed 



u 



Onsager 
|2/3=log(l+v^) 



-V2 = Ml(/?e 



(22) 



where u^'^'^s^'' {(3) = ~ctgh{2(3) 



1 + 



K 



2tanh(2/5)^ - 1, + 



J2 



K'" = 1, Ki is the elliptic function of the first kind [|l|], and ui{(3) = — tanh {(3)- 
1. For d > 2 Eq.(pTD remains unproven similarly to Eq.(^. 

To conclude, there are many unanswered questions and much more work 
is to be done, but we feel that it is certainly worth to undertake this effort. 
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calculating low moments of the transfer matrix. I would like to thank him 
for numerous and stimulating discussions. I also thank E. Seller and A. Sokal 
for instructive discussions. 
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Figure Captions 

Fig.l. Construction of the transfer matrices T and T2 for the two dimen- 
sional Ising modeL Periodic boundary conditions are understood. 
Fig.2. Same as Fig.2 but for the triangular lattice. Bonds implied by the 
periodic boundary conditions are shown explicitly. 
Fig. 3. Same as Fig.3 but for the two layer Ising model. 
Fig.4. Characteristic function for the two layer Ising model. 
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Fig. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9312025vl 



